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Abstract 

We show that the condition for the appearance of quantum chaos 
(Wigner-Dyson distribution of energy eigenvalues, gaussian-random 
energy eigenfunctions) in a dilute gas of many hard spheres is A <C £, 
where A is the wavelength of a typical particle and i is the mean free 
path between collisions. For fermions with Fermi wavelength Xp £, 
this implies that all energy eigenstates, including the ground state, are 
chaotic. Physical implications are discussed. 
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Quantum systems which are classically chaotic |l| are known to exhibit several interesting 
and universal features, such as Wigner-Dyson statistics for energy eigenvalues |2j , gaussian- 
random energy eigenfunctions || , and quantum expectation values which approach classical 
microcanonical averages . These features can then be used to understand certain physical 
phenomena, such as the statistics of fluctuations in the conductance of quantum dots 
0. We will refer to these features collectively as "quantum chaos". Quantum chaos is 
expected to appear only for sufficiently high energy eigenvalues; an important problem is 
the determination of the energy scale which marks the (possibly gradual) onset of quantum 
chaos. Here we consider this problem for a dilute gas of hard spheres, a many-body system 
which is (almost certainly 0) classically chaotic. 

In general, the presence of quantum chaos at an energy E requires two distinct conditions 
to be satisfied. The first is that 

A</ry s i ow , (1) 

where A is the mean spacing of energy eigenvalues near energy E, and 7 s i ow is the rate of 
the slowest classical process which is physically relevant. For a system with n > 2 degrees 
of freedom, condition ([]]) is satisfied in the formal limit of % — > with E fixed, since in that 
limit A otTi n . If, on the other hand, we consider H to be fixed, then ([]]) implies a lower bound 
on the energy where we can expect the signatures of quantum chaos to appear. Condition 
(1) can be understood heuristically by considering the corresponding time scales (e.g., @). 
Any relevant classical process must occur before the Heisenberg time th = 2irh/A] this 
is because th is the time at which all observables begin to oscillate quasi-periodically, an 
inherently quantum mechanical phenomenon which shuts off any later classical effects. 

For example, for a single particle in a random potential in d dimensions (representing, 
e.g., an electron in a disordered metal), 7 s i ow ~ D/L 2 is the Thouless rate [fj; here D ~ £v 
is the classical diffusion constant, £ is the mean free path, v = (2E/m) l l 2 is the speed of the 
particle, and L ^> i is the linear size of the system. For the Sinai billiard in d dimensions, 
in which a particle moves ballistically in a cubic box with a single spherical scatterer in the 
center, 7 s i ow ~ v/£ ~ a v/L , where v is the speed of the particle, L is the linear size of 
the box, a is the radius of the spherical scatterer, and I ~ L d /a d ~ l is the mean free path 
between collisions of the particle with it; here we have assumed o C i, and hence L <C t. 
In both examples, A ~ (X/L) d E, where A = 2irh/mv is the wavelength of the particle. 

For the example of a particle in a random potential, condition (1) leads to A d_1 <C £L d ~ 2 . 
When d — 1, this condition is independent of A, and then it is not satisfied (since, by 
assumption, i <C L in this case). For a particle in a random potential, quantum chaos 
corresponds to the existence of extended (rather than localized) energy eigenstates [§], so 
(|I|) correctly predicts that there are no extended states when d = 1. It will turn out that 
the second condition, specified by Eq. (||D below, is more stringent than ([!]) when d > 2, and 
so we defer further discussion of this example. 

For the example of the Sinai billiard, (|l|) leads to A <C a. This condition is also necessary 
in order to have large phase shifts when a particle with wavelength A is scattered from a 
sphere of radius a in free space. The appearance of quantum chaos only for A <C a is also 
supported by some numerical computations |J. 

The second condition which must be satisfied for the appearance of quantum chaos is 

&7fest < E , (2) 
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where E is the total energy of the system and 7f ast is the rate of the fastest classical process 
which is physically relevant. Like (|l|), this condition is satisfied in the formal limit of h — > 
with E fixed, and like (|I|) it can be understood heuristically by considering the corresponding 
time scales. For a quantum state whose energy uncertainty is no larger than the energy itself, 
h/E represents the minimal time uncertainty in any physical process. A classically relevant 
process which takes less time than h/E will not be reflected in the quantum theory. 

In both of our examples, 7f ast ~ v/£ is the inverse of the mean free time; condition (^|) 
then implies A <C i. In the example of the Sinai billiard, this is weaker than the condition 
A < a which follows from ([TJ). In the example of a particle in a random potential, A <C I is 
well-known from the theory of weak localization as a necessary condition for the existence 
of extended states when d > 3 (e.g., fLOp . The critical case of d = 2 is more complicated: 
all energy eigenstates are localized, but the localization length is given by £ ~ exp(c£/A)£, 
where c is a numerical constant. If £ ^> L, then the states are effectively extended and the 
signatures of quantum chaos will be present. 

The most rigorous derivation of either (p]) or (|2|) follows from the recent work of Andreev, 
Agam, Simons, and Altshuler (AASA) [|11|], who computed statistical properties of energy 
levels from first principles, using only generic classical features of chaotic systems, and 
making reasonable (but not strictly proven) assumptions about the form of the classical 
limit of certain functional integrals. It is clear that the approach of AASA extends in a 
straightforward way to properties of eigenfunctions and matrix elements. For a classical 
system with a Perron- Frobenius spectrum (that is, the eigenvalues of the classical time 
evolution operator when it is restricted to distributions on phase space that are suitably 
well-behaved at all times) that has a gap between the lowest eigenvalue of zero and the first 
nonzero eigenvalue 71, AASA find that 

7siow = Re7i. (3) 

For a system with no gap, 71 can be identified with the inverse of the time T\ at which classical 
correlation functions begin to decay as a power of the time (instead of exponentially). The 
heuristic estimates of 7 s i ow which we made in our two examples are consistent with Eq. @. 

The value of 7f ast in the AASA formalism follows from the requirement that the Weyl 
symbol of a product of two operators is equal (approximately) to the product of their 
individual Weyl symbols. In particular, this must apply to products of the hamiltonian 
with itself. Given two operators A and B with Weyl symbols Ay/(x) and Byy(x), where 
x = (q,p) are canonical coordinates on phase space and 

A w {x) = J d d q'e-^'/ h (q + \q'\A\q - \q') , (4) 

the Weyl symbol of their product AB can be expressed as [ jl^l 

(AB) w (x) = exp [(in/2)Vr J-V a ] A w (x 1 )B w (x 2 )\ Xl = X2=x , (5) 

where J = ((0, —I), (1, 0)) is the simplectic matrix. To find 7f ast , we must compute the 
discrepancy between (H 2 )w(x) and [Hw(x)} 2 , and demand that it be small. We begin 
by recalling that, in a chaotic system, unstable periodic orbits are dense in phase space, 
and that the number of them with period T is proportional to exp(/i top T), where h top is the 
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topological entropy. Given an unstable periodic orbit, we can choose phase space coordinates 
x = (E,t,X) where E is the energy, t measures length along the orbit, and X represents 
coordinates in the Poincare section of the orbit. For small X we can write the Weyl symbol 
of the hamiltonian as 

H w (x) =E + \X-XX, (6) 

where A is a diagonal matrix whose entries are the Lyapunov exponents Aj, which are real 
and come in pairs with equal magnitudes and opposite signs. Applying (|5[), we find 

{H 2 ) w {x) = [H w {x)f - \h 2 Ei A, 2 + 0{h A ) . (7) 

Thus we must have Ti 2 J2i A 2 <C E 2 in order for the symbolic manipulations of AASA to be 
valid. 

A similar but somewhat stronger condition follows from a closely related point of the 
AASA analysis; given any function on phase space fw( x ) which is smooth on the scale set 
by h, AASA assume that (Hf) w (x) is approximately equal to H w (x)fw(x). Let us consider 
the specific case of 

f w (x) = l + (a/h)X-K-X + ... , (8) 

where K = ((/, 0), (0, — /)), and a is a numerical constant (which must be small if fw{x) is 
to be slowly varying). We also write A = ((A + , 0), (0, — A + )), where A + is a diagonal matrix 
whose entries are the positive Lyapunov exponents A^~. Then we find 

(Hf) w (x) = H w (x)f w (x) + \ahY, t A+ + 0(h 2 ) . (9) 

Thus we must have hJ2i^t ^ E- The sum of the positive Lyapunov exponents is the 
Kolmogorov- Sinai entropy hxs, and so we can make the identification 

7fast = h KS . (10) 

Of course, by choosing other forms for the matrix K, we can equally well conclude that we 
must have hJ2i a i^i ^ E for any coefficients «j which are each of order unity or smaller. 
The Kolmogorov- Sinai entropy /iks is simply the largest linear combination of this type. 

Eq. (p~0| ) also follows from an analysis of the validity of the "diagonal approximation" 
13| to the semiclassical periodic orbit expansion of, e.g., the spectral density correlator. 



When valid, the diagonal approximation can be used to derive some of the signatures of 



quantum chaos ||14|| . Without employing any resummation techniques, the periodic orbit 
expansion converges only if the energy is given an imaginary part whose magnitude exceeds 
fr(h t0 p — |/iKs) jl5f . When E is kept strictly real, the diagonal approximation relies on large 



E to produce rapidly varying phases which then yield large cancellations. These can be 
expected to occur only if E 3> h(h top — ^h KS ). Since in general h top > /iks, we once again 
find the condition hh^s <C E. 

We now turn to our main subject, the application of (|I|) and (Q) to a gas of ^> 1 hard 
spheres, each with radius a, in a box with edge length L in d dimensions. We will take the 
gas to be dilute, so that Na d L d . The length of the mean free path between collisions of 
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one particular sphere with any other is then i ~ L d / N a d ~ l , which is much larger than a; we 
will, however, assume that f < L. 

We must now determine, as functions of the total energy E, the key quantities 7 s iow, 
7f ast , and A. The low lying modes of the Perron- Frobenius operator are associated with 
diffusive processes, and so we expect 7 s i ow ~ D/L 2 ~ £v/L 2 , where v = (2mE/N) 1 ^ 2 is the 
typical speed of a single particle. Thus, in the thermodynamic limit of iV — > oo with E/N 
and L d /N fixed, 7si ow goes to zero like N~ 2 / d . However, in a many-body system, the level 
spacing A is related to the thermodynamic entropy S via A ~ e~ s E. Since entropy is an 
extensive quantity, A is exponentially small in N as N — > oo. Therefore condition (|T|) is 
always satisfied for large N, independent of E. 

To evaluate condition ([|), we need the Kolmogorov- Sinai entropy of the hard-sphere gas; 
it is given by 



h KS ~ (Nv/£) ln(£/a) , (11) 

where again v = (2mE / N) 1 ^ 2 is the typical speed of a single particle. If define the single- 
particle wavelength A = 2irh/mv, we then find that condition (g) yields 

a<< rW (12) 

This, then, is the requirement for the appearance of quantum chaos in a dilute gas of many 
hard spheres. 

Let us eliminate a possible loophole in the preceding argument. We have assumed that 
the spectrum of single-particle wavelengths contained in an energy eigenstate is characterized 
by a single length scale A ~ (h 2 N /mE) 1 / 2 . Fortunately, in the energy range where the 
signatures of quantum chaos appear, this is correct: the distribution of single-particle 
wavelengths in a chaotic energy eigenstate is thermal, with the characterizing temperature 
related to the energy eigenvalue by the usual thermodynamic relation between temperature 
and energy [ I8| . Of course, which thermal distribution we get (Maxwell-Boltzmann, Bose- 



Einstein, Fermi-Dirac) depends on the assumed statistics of the particles, an issue we have 
ignored up to now. 

Here we will consider the case of fermions, which, for sake of definiteness, we will assume 
to have spin one-half. (Spin degeneracy will, however, affect only various numerical factors 
which we will not evaluate explicitly.) The Weyl symbol calculus can be modified to reflect 
Fermi statistics |]L9|| , as can the periodic orbit expansion flTEIJZUfl . Therefore all of our previous 



analysis is still applicable. The ground state of the system is characterized by a Fermi 
wavelength A^ ~ L/N 1 ^; note that, for a dilute gas, a < Ai? <C Since \p <C £, 
the ground state (and, indeed, every energy eigenstate) apparently meets our test for the 
appearance of quantum chaos. 

Our claim of quantum chaos in the ground state must be reconciled with the standard 



theory of the normal Fermi liquid [^TJ , in which quasiparticle excitations have weak couplings 
characterized by the small parameter a/Ai?. The key point is that quasiparticles do not 
represent energy eigenstates; this is clear from the facts that quasiparticle energies have an 
imaginary part, and that their number is not conserved. Instead, quasiparticles represent 
superpositions of energy eigenstates with particularly simple behavior under time evolution. 
A quasiparticle is analogous to a wave packet in the Sinai billiard, narrow in both position 
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and momentum, constructed by appropriate superposition of high energy eigenfunctions. 
Such a packet will move simply (that is, ballistically) for a long time. This is true even 
though the individual energy eigenfunctions of which it is composed are themselves chaotic. 
The difference from the case of the normal Fermi liquid is that there we do not require high 
energies for the occurrence of quantum chaos. 

Another point is that the results of diagrammatic perturbation theory (for the ground 
state expectation values of simple operators) imply that the exact ground-state wave function 
can be written as an asymptotic series of the form 

Mi) = 4°\q) + (a/x F )4 1] (<i) + -.., (13) 

where ifj^ (q) is the usual Slater determinant of eigenfunctions of a single particle in a box. 
Quantum chaos implies (among other things) that the distribution of the values of ipo(q) 
over the configuration space is gaussian 0,^,17]. This can be reconciled with ([13]) only if 



the same is true of the zeroth-order approximation ip^iq). A gaussian distribution of the 
values of ip^iq) implies, and is implied by, 



L J 10 for k odd 

where we have set L = 1 for notational simplicity. Fortunately, Eq. (0) can be shown to be 
valid for large N. This is more easily demonstrated if we put periodic boundary conditions 
on the box (instead of hard walls) so that the single-particle eigenfunctions are plane waves. 
Then, if we keep only those terms on the left-hand side of ([14]) in which each single-particle 
eigenfunction Xp( x ) = e ip ' x / s appears together with x_ p (x), the result is the right-hand side 
of (|T4|). The remaining terms on the left-hand side are much more numerous, but most of 
them integrate to zero, and those that do not can have either sign. Assuming that these 
signs can be treated as statistically random, one can show that the net contribution of the 
remaining terms to the right-hand side of (0) is exponentially small in N. 

The appearance of quantum chaos at all energies in the normal Fermi liquid does have 
physical consequences for quasiparticle behavior on long time scales. It predicts that thermal 
equilibrium will eventually be attained, even in the complete absence of any coupling to an 
external environment, and beginning from an arbitrary initial state [pTS| . Thus quantum 



A' 



(k — 1)!! for k even 



chaos can serve as the fundamental dynamical underpinning of the usual phenomenological 
transport equations (which also predict, in more detail, an approach to thermal equilibrium). 
In this regard, it is interesting to note that the condition A <C I also arises during a careful 
derivation of a Markovian transport equation for a dilute gas with a smooth, short-range 
interaction potential between particles ||22|| . (This formalism, however, is applicable only to 



a highly restricted class of initial states.) 

As another possible physical consequence, let us note that a dilute gas of hard-sphere 
fermions in a box is a poor model of a nucleus. Nevertheless, if we extrapolate the present 
result to a dense gas with a weak long range attraction (replacing the box) but a short range 
repulsion, it may be that quantum chaos also extends down to the ground state in large 
nuclei . 



To conclude, we have explored the conditions for the appearance of quantum chaos in a 
dilute gas of hard spheres, and found that they imply that a typical particle's wavelength A 
must be much less that the classical mean free path £. This condition is satisfied by all energy 
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eigenstates, including the ground state, if the particles are fermions. This result is consistent 
with the standard theory of the normal Fermi liquid. It can help to put phenomenological 
transport equations on a firmer theoretical foundation, and may have implications for the 
properties of low-lying energy eigenstates in large nuclei. 
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